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The transition from a spectrum filter (resonant transmission) to a polarizer (broadband transmission) for 
TM polarized light is observed in a metallic nano-slit array as period is decreased. A theoretical model is 
developed and shows that the spectrum filter behavior is caused by the coupled slit/grating resonance. With 
decreasing period, the slit resonance is decoupled from the grating resonance, which then dominates the 
transmission spectrum and broadens the transmission peak. With further reducing period, the slit 
resonance diminishes and the peak spectrum transforms to a broadband transmission. This effect is the 
basis for the operation of wire grid polarizers. The transition is explained by the change of the impedance to 
the incoming wave. 

Sub -wavelength metallic nano-slit arrays have attracted much attention due to the extraordinary resonant 
transmission of transverse magnetic (TM) waves 1 " 12 . Based on the resonance behavior, metal slit arrays have 
been studied as polarizing spectrum filters 13,14 and utilized in band-pass photodetectors 15 . Recently they 
become a promising candidate of the structural color filters in flat-panel displays and colored imaging sensors 1617 . 
On the other hand, of a similar structure, metal wire grids have been extensively utilized as polarizers, with the 
characteristic high TM transmittance over broadband wavelength range 18 " 23 . A question naturally emerges: why 
does the similar kind of structures (metal slit arrays) have drastically different transmission characteristics? This 
work aims to reveal the mechanisms that are responsible for both effects. Especially we present experimental 
results and analysis to show the transition from a spectrum filter (resonant narrowband) to a polarizer (flat 
broadband) for TM waves in metal slit arrays. 

A general observation is that the metal slit array based polarizers (or more commonly referred to as wire grid 
polarizer, WGP) usually have a smaller period than the spectrum filters. For example the typical period of a WGP 
working in the visible band is less than 150 nm 20 . Therefore, we study the evolution of the transmission spectrum 
of a metallic slit array as a function of its period. Experimentally we observed transition from a spectrum filter to a 
polarizer with decreasing period. In general there are two types of resonances associated with the structure: 1) the 
grating resonance, which is due to constructive interference of the scattered light by the periodic grating struc- 
ture 24 " 27 and 2) the slit resonance from a localized Fabry-Perot resonance of the slit waveguide mode 28 " 32 . Many 
effects involve the interplay of both effects. Our theoretical analysis shows that the transition from resonant 
transmission in a spectrum filter to broadband transmission in a polarizer involves two processes: 1) decoupling 
of the slit resonance from the grating resonance and 2) diminishing of the slit resonance due to cancellation of 
induced polarization charges at the slit entrance. The slit resonance itself has a broad profile, which is confirmed 
by the transmission measurement of a slit array with very large period. Through the analogy between a metal slit 
array and a dielectric slab (following Ref. 33), an effective refractive index was derived which describes the 
diminishing process of the slit resonance, eventually giving rise to the broadband transmission of polarized 
TM light as in a typical WGP. In addition, we studied the angular dependence of the device as it changes from 
a spectrum filter to a polarizer and we found that the device becomes less angular dependent with decreasing 
period during the transition process. 

Results 

Decoupling of the coupled slit/grating resonance with decreasing period. Metallic nano-slit arrays are 
experimentally realized by creating perforating slits in an Al film on a fused silica substrate through standard 
e-beam lithography and reactive ion etching. Four periods (d = 200 nm, 300 nm, 400 nm and 500 nm) with the 
same slit width (a ~ 110 nm) and the same slit height (/z — 170 nm) are fabricated. The sketches and the SEM 
images of two samples with a period of 400 nm and a period of 200 nm are shown in Fig. 1. Firstly, the samples 
with four different periods are all polarization dependent. Typically both spectrum filters and polarizers require 
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Figure 1 | Sketches and the corresponding SEM images of two metallic nano-slit arrays with different periods, (a) and (c) show a fused silica supported 
Aluminum slit array with a larger period (400 nm) serving as a spectrum filter, (b) and (d) show the similar structure with a smaller period (200 nm) 
serving as a polarizer. Both these two structures have the same slit width of —110 nm, and the same slit height of — 170 nm. 



sub -wavelength scale metallic gratings, the incident TE wave is 
beyond the cut-off wavelength of the waveguide formed between 
the metal slit walls, therefore the TE waves are mostly reflected. So 
we focus on the behavior of TM polarized light. The transmission 
and reflection spectra together with photos of the four samples 
illuminated by a TM wave at normal incidence are shown in Fig. 2 
(b) and (c), respectively. At a large period (e.g. 500 nm or 400 nm), 
the resonant spectrum is prominent, giving rise to the distinctive 
transmission and complementary reflection colors. Resonant trans- 
mission peak or reflection dip is observed with asymmetric line 
shape. With decreasing period, the resonant peak or dip broadens 
and also shifts to shorter wavelengths. At the period of 200 nm, the 
peak is so broad that the spectrum looks almost flat for the whole 
measurement range. Accordingly, a white transmission and a dark 
reflection for TM light are shown in the photos. In this experiment, 
the transition from a spectrum filter behavior (i.e. a resonant 
narrowband transmission) to a polarizer behavior (i.e. a broad- 
band flat transmission) is clearly realized by simply varying the 
period of the metallic grating. 

We used an analytical model based on mode expansion (see the 
Supplementary information online) to understand the observed 
effects. Besides the basic parameters of the slit array structure (such 
as slit width, slit height and period), the incident angle, the refractive 
index of the substrate and that of the slit filler can all be tuned in this 
model so that different measurements can be compared with theory 
by tuning the corresponding parameters in the model. Although the 
model assumes the metal to be a perfect conductor as an approxi- 
mation, it captures most of the experimental features and provides a 
reasonable physical picture. Similar models have been established to 
describe the optical response of a grating 1 ' 2 ' 5 ' 25 27 ' 29 31,34 36 . We will 
discuss the effect of using real metal in the experiment later in this 
paper. 

Due to the sub -wavelength scale of the slit array, the measured 
transmission is of the zeroth order (see the Methods section). At 
normal incidence and with air slits, the zeroth order transmission 
is derived by the model as 



where is the dielectric constant of the substrate; a denotes the slit 
width; d denotes the period; a + and a~ are the coefficients of the 
fundamental slit waveguide modes propagating forwards and 



backwards, respectively. Following the experimental conditions, 8^ 
is chosen to be 2.10 for a fused silica substrate; a and d are chosen to 
be 1 10 nm and 170 nm, respectively. A 2D contour plot of the zeroth 
order transmission T 0 versus wavelength and period is shown in 
Fig. 2 (a). The red zone in the color map corresponds to the resonant 
transmission peak. The oblique dark lines (marked out by dashed 
white lines) represent the different diffraction orders of the grating 
resonances. They are divided into two sets due to the higher refractive 
index of the substrate than air. The two lines marked by N' = 1 and 
N' = 2 represent the first and the second order grating resonance at 
the grating/substrate interface while N = land N = 2 represent the 
diffractions at the grating/air interface. The pronounced transmis- 
sion resonance in the 2D contour plot is due to the coupling of the 
grating resonance with the slit resonance 1 ' 3 " 6 ' 8 ' 9 . The slit resonance is a 
result of the standing wave formation due to the two counter-prop- 
agating waves represented by a + and a~, i.e. a Fabry-Perot resonance. 
Such slit resonance is independent of the period due to its localized 
characteristic. In addition, the slit resonance (marked by the band 
between the two vertical dash lines in Fig. 2 (a)) is much broader than 
the grating resonance (due to the low impedance contrast, as will be 
discussed further later in the paper). When the grating resonance and 
the slit resonance overlap in spectrum at certain specific periods, the 
two resonance modes couple to each other and produce the familiar 
anti-crossing behavior. The resonant peak which is red-shifted with 
increasing period is actually the lower branch of the anti- crossing. 
The higher branch is coupled to a higher order grating resonance so 
that it is not observable due to the finite numerical aperture used in 
our experiment. The interference between a broad resonance (the slit 
resonance) and a sharp resonance (the grating resonance) gives rise 
to the typical Fano-type line shape, as confirmed by the asymmetric 
shape of the resonant peak in the measured transmission or the dip in 
reflection. The transmission vs. wavelength presented in this 2D 
color map for the period used in the experiment agrees with the 
measurement data (see the Discussion section and the Supplemen- 
tary information online). At the period of 500 nm, the grating res- 
onance overlaps the central wavelength of the slit resonance so that a 
strong coupling is formed. Due to the high-Q characteristics of the 
grating resonance 3 ' 4 , the coupled slit/grating resonance is much 
sharper than the intrinsic slit resonance, making it attractive to the 
spectrum filter application. With decreasing period, the grating res- 
onance gradually moves away from the slit resonance towards 
shorter wavelengths. Consequently, the slit and grating resonance 
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Figure 2 | Transmission and reflection measurements on the slit arrays with four different periods, in comparison with a period dependent 
transmission spectra based on the theoretical model, (a) shows the 2D plot of a period dependent zeroth order transmission spectra based on the 
theoretical model. The oblique dashed lines mark the grating resonances of different orders. The band between the two vertical dashed lines marks out the 
slit resonance, (b) and (c) show the measured transmission and reflection spectra (black dots) together with RCWA simulations (red curves) of the fused 
silica supported Aluminum slit arrays at four different periods (200 nm, 300 nm, 400 nm and 500 nm). Slit width ( — 110 nm) and slit height 
( — 170 nm) remain invariant for different periods. Insets show the transmission and the reflection photographs of the four samples. Each square is 
200 urn X 200 urn. 



become separated, resulting in a broad transmission peak due to the 
slit resonance and a sharp transmission dip due to grating resonance. 
Decoupling of the two resonances with decreasing period is also 
manifested by the reduced angle dependence as shown in Fig. 3. As 
well known, the grating resonance is strongly dependent on the 
incident angle. As shown in Fig. 3, with an increasing incident angle, 
each grating resonance splits into two branches due to the phase 
matching requirement. The branch shifting to shorter wavelengths 
represents the forward diffracted light whose k x has the same sign as 
that of the incident light. The branch shifting to longer wavelengths 
represents the backward diffracted light whose k x has the opposite 
sign to that of the incident light. On the other hand, the slit resonance 
is intrinsically independent of angle 17 . At the period of 400 nm, the 
backward branch from the first order grating resonance at the grat- 
ing/substrate interface (marked as b') meets the slit resonance, form- 
ing a strong coupling between the two. The hybrid mode takes on the 
characteristics of both, e.g. the overall transmission peak becomes 
angle-dependent as in the case of grating resonance. The theoretical 
plot with an extended angle range (0-20 deg) (Fig. 3 (c)) shows that 
the transmission peak shifts to longer wavelengths with increasing 
angle. When the period decreases to 300 nm, the b' diffraction 
branch is farther away from the slit resonance so that the two reso- 
nances decouple and consequently the transmission peak due to the 
slit resonance becomes less angle-dependent. At the period of 
300 nm, an angle-insensitive range from 0 to 8 deg is observed in 
the theoretical plot in Fig. 3 (d). The angular dependence behavior 
will be further discussed in the last part of the Results section. 

The intrinsic slit resonance. From the period dependent spectra 
(Fig. 2) and the angle dependent spectra (Fig. 3), the slit resonance 
is almost decoupled from the grating resonance at the periods below 
300 nm. In this case, the transmission peak is mostly dominated by 
the slit resonance. In order to reveal the intrinsic slit resonance, we 
fabricated the metal slit arrays with very large periods to emulate the 



case of a single-slit. The measured transmission at such large periods 
(e.g. 5 |im for Fig. 4 (c)) is very weak, however a broad peak in 
transmission accompanied by a broad dip in reflection are still 
discernable (Fig. 4 (d)). The oscillatory data superimposed on the 
broad background is well above the noise level in the measurement. 
Although the zeroth order transmission is experimentally accessible, 
it is not a good representation of the slit resonance, especially at a 
large period because many diffraction orders of light transmitting 
through the slit cannot be collected by the finite aperture in the 
experiment. Therefore, we introduce the slit power enhancement 
factor, which is the ratio of the power flow inside slits ((S s i it )) to 
that of the incident light ((Si)). At normal incidence, it can be 
expressed as 



The light power inside the slits is directly influenced by slit 
resonances so that the slit power enhancement factor is more 
sensitive to slit dimension itself than to grating structure. In 
experiment, even though the light power inside the slits spreads 
into different diffraction orders, the slit power enhancement factor 
encompass all the diffraction orders and therefore remains high at 
large periods. By plotting this quantity versus wavelength and period 
in a 2D contour map, the slit resonances are clearly revealed against 
the grating resonances at all periods, as shown in Fig. 4 (b). The two 
vertical stripes in the 2D contour map of (S s i it )/(S t ) (Fig. 4 (b)) are due 
to the first and the second order Fabry-Perot resonance in a single 
slit, whose bandwidths can be read out from Fig. 4 (a) as —450 nm 
and —100 nm, respectively, and they are obviously independent of 
period; whereas the dense array of oblique lines are due to the grating 
resonances, and are linearly dependent on the period as expected. 
The slit resonance does cut through a number of grating resonances 
as shown in Fig. 4(b), forming the coupled hybrid modes as discussed 
above, which are responsible for the small oscillations superposed 
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Figure 3 | Angle- dependent transmission measurement at two periods together with the theoretical results, (a) and (b) show the measured angle- 
dependent transmission spectra for fused silica supported Aluminum slit arrays with a period of 400 nm and a period of 300 nm, respectively. The slit 
width of the structure is —110 nm and the slit height is —170 nm. (c) and (d) show the angle-dependent transmission spectra by the theoretical model in 
an extend angle range. The branches marked out by the thin dashed lines are due to grating resonances, /and b mark the forward and the backward 
branches, respectively, of the grating resonance at the grating/ air interface; / and b' the forward and the backward branches at the grating/substrate 
interface. The thick dashed lines marks the transmission peak. 



over the broad slit resonance in the measured spectrum (black 
squares in Fig. 4 (c) and (d)), and qualitatively agree with the 
Rigorous Coupled Wave Analysis (RCWA) simulation results (red 
curves in Fig. 4(c) and 4(d), also see the Discussion section). 

Dimished slit resonance with further decreasing period— tran- 
sition to broadband polarizer. An illuminated metal slit with a 
deep sub-wavelength dimension can be considered as a dipole 37 39 . 
When driven by a TM wave with the direction of the electric field 
crossing the slit, polarization charges are accumulated at the corners 
and consequently distort the nearby field. At either the entrance or 
the exit side of the slits the incident or the outgoing plane waves 
should couple into or out of the waveguide modes in the slits, 
which creates an impedance mismatch at the two interfaces. The 
impedance difference leads to partial light reflection from the two 
surfaces. Such impedance mismatch and the constructive field 
buildup of slit waveguide modes establishes a Fabry-Perot type 
resonance— this is the nature of the slit resonance discussed above. 



In the case of an array of slits, the separation between each slit plays 
an important role in deciding the impedance mismatch and 
consequently the light confinement. For a 400 nm period, as 
shown in Fig. 5 (a) p400, the induced dipoles are separated from 
each other by sufficient distance. Thus, the light field around each 
dipole is distorted, leading to a large impedance mismatch between 
the incident plane wave and the slit waveguide mode. Thus, light is 
strongly reflected at both ends of the slit, which is desirable for 
forming high-Q F-P resonance modes. When the period decreases 
to 200 nm while the slit width remains 110 nm, the induced dipoles 
are adjacent to each other end to end as shown in Fig. 5 (a) p200 so 
that the accumulated charges tend to cancel out. Thus, the light field 
is only slightly distorted. Therefore the impedance mismatch to air is 
significantly reduced, leading to a much weaker F-P resonance and 
significantly broadened transmission "peak". Indeed the broadband 
polarizer is achieved when the slit resonance becomes very weak so 
that high transmission is maintained across a wide wavelength range 
of the slit resonance. 
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Figure 4 | Identifying the intrinsic slit resonance, (a) Spectrum plot of the slit power enhancement factor ({S s i it )/{S t )) at the period of 1 mm based on the 
theoretical model. Inset illustrates the definition of the slit power enhancement factor as a ratio of the power flow inside the slit to that of the incident light, 
(b) Period-dependent spectra of the slit power enhancement factor. The white dashed line marks out the spectrum at the period of 5 um to compare with 
the measured data, (c) and (d) Measured transmission and reflection spectra (black dots) together with RCWA simulations (red lines) of the fused silica 
supported Aluminum slit array with a period of 5 (am. An eye-guide line (yellow dashed) is plotted in (c) and (d) to reveal the broad slit resonance. The slit 
width is — 110 nm and the slit height is — 170 nm. The reflection spectrum is in agreement with the simulation while the transmission is about three times 
higher. It is probably due to the high magnification of the microscope objective (40 X ) that was used in transmission measurement in order to enhance the 
signal by collecting several diffraction orders instead of only the zeroth order. Insets show the transmission and the reflection photographs of the slit array. 



Finally we provide an analytical approach to describe the above 
mechanism. The metal film with perforating slit array can be mapped 
into a dielectric slab with an effective refractive index (n e ff) and an 
effective thickness (h e jfi by comparing either the reflection coefficient 
or the transmission coefficient of a free-standing metal film with that 
of the dielectric slab. In this sense, the slit resonance of the metal film 
with perforating slit array corresponds to the F-P resonance of the 
effective dielectric slab. This method follows that prescribed in Ref. 
33. Assuming the metal to be a perfect conductor, n e ff and h e jj are 
expressed as 

1 



+ oo 

Ea ( a\ 
-sinc^-j 



1 



(3) 



and h e ff = hln e ff. Through this analogy, the impedance mismatch is 
represented by n e ff, which according to equation (3) depends on two 
ratios, the ratio of slit width to period (aid) and the ratio of wave- 
length to period (XI d). In the long wavelength limit (i.e. kld^> oo), the 
higher orders in the sum of equation (3) only introduce a small 



imaginary correction, and then n e ff approximately equals to dla, 
which is proportional to period. Consequently, with decreasing per- 
iod the reflectance reduces as shown in Fig. 5 (d). It is well known that 
the reflectance of a film determines the finesse of the F-P etalon made 
of this film. With reducing reflectance, the F-P resonance of the 
effective dielectric slab (i.e. the slit resonance of a metal slit array) 
diminishes so that the transmission become broadband and flat. At 
smaller wavelengths, although n e jj becomes a complex number as 
shown by two examples (600 nm and 1000 nm) in Fig. 5 (c), the 
corresponding reflectance at the metal slit array still decreases with 
decreasing period (Fig. 5 (d)), which confirms the diminishing pro- 
cess of the slit resonance. 

Angular dependence. As the metal slit arrays transition from a 
spectrum filter to a polarizer due to decoupling of the slit/grating 
resonance, the angular dependence reduces simultaneously. Thus, a 
polarizer made of the metal slit array structure becomes increasingly 
angle-independent with reducing periods. On the other hand, an 
angle-insensitive spectrum filter is desirable in many applications. 
As discussed, the angular dependence can be reduced by decreasing 
the period. However the transmission peak broadens significantly 
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Figure 5 | Illustration of the slit resonance diminishing with period decreasing, (a) Simulated x- component E field distributions and electric field lines 
of two metal slit arrays at their own resonant wavelengths, i.e. the period of 400 nm at its resonant wavelength of 687 nm; the period of 200 nm at its 
resonant wavelength of 507 nm. As revealed by the calculated electric lines, the light field near the slit openings in the 400 nm-period array is distorted 
more severely than that in the 200 nm-period array. Accordingly, the calculated x-component of electric field (E x ) reveals a stronger on-resonance light 
field inside the slits of the 400 nm period than that of the 200 nm period. The max value of the colorbar is 2.44 X 10 5 V/m; the min value is —2.90 X 
10 5 V/m. (b) Sketch of the analogy between a metal slit array layer and a dielectric slab, (c) Evolution of the effective refractive index {n e fj) with period at 
three different wavelengths (i.e. 600 nm, 1000 nm and infinite) plotted in a complex plane with the horizontal and vertical axis to be the real and 
imaginary part of n e ff. Period increases along each curve in the direction of the arrows. Three periods, i.e. 200 nm, 300 nm and 400 nm are marked out by 
squares, circles and triangles, (d) Evolution of the reflectance from the metal slit array with period. The accurate reflectance (R0_acc) is in squares and the 
reflectance calculated through the effective index (R0_eff) is in solid lines. At an infinitely large wavelength, these two quantities are identical. The 
expression of R0_acc can be found in the Supplementary information online. 



due to two reasons: 1) the decoupling of the sharp grating resonance 
from the broad intrinsic slit resonance increases the band width of 
the transmission peak. 2) Reduction in the impedance mismatch 
between the light mode inside and outside the slits further broa- 
dens the slit resonance. Therefore, in order to design an angular 
insensitive spectrum filter using this structure, the bandwidth of 
the slit resonance should be reduced. This can be done by either 
shrinking the slit width or filling the slit with a high index 
dielectric, so that the Q factor of the slit cavity is enhanced and 
resonance bandwidth reduced. 

Discussion 

Our theoretical model assumes the metal to be a perfect conductor, 
which is an approximation of the real metal used in the experiments. 
The validity of this approximation in describing the transition from a 
spectrum filter to a polarizer in a real metal slit array is studied. As a 
direct comparison, the transmission spectra predicted by the theor- 
etical model, the measurement results and the simulation by using 
RCWA based on real metal index are plotted together in Fig. 6 (a). 
The period of the metal slit array is 400 nm. The slit width and the slit 



height are 110 nm and 170 nm, respectively. Clearly the theoretical 
model based on ideal metal and RCWA simulation based on real 
metal both reproduce the main features observed in the experiment, 
i.e. the grating resonances (represented by the two dips) and the slit 
resonance (represented by the broad peak). There are also discrep- 
ancies in the transmission amplitudes and resonance peak positions: 
1) the spectrum of the Al slit array (either the measurement or the 
RCWA simulation) shows a lower transmission peak at the slit res- 
onance than the theoretical prediction based on ideal metal, which is 
due to the optical loss in Al; 2) both the slit resonance and the grating 
resonances are red- shifted in the spectrum of the Al slit array with 
respect to the theoretical model. This can be explained by the differ- 
ent types of propagating waves for the real metal and the ideal metal 
case, in particular the difference in propagation constants. Simply 
put it, the grating resonance (Fig. 6 (b)) in the ideal metal case is due 
to the diffracted wave along the surface direction having free space 
momentum, while for real metal it is due to the coupling to the 
surface plasmon polariton (SPP) wave, which is known to have 
higher momentum due to its special dispersion, and therefore 
occurs at longer wavelength (see Supplementary Fig. S3 (c) online). 
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Figure 6 | Comparison between a slit array made of a perfect metal and that made of a real metal, (a) Transmittance of a slit array made from a perfect 
metal (theoretical prediction: blue curve) or Al (measured data: black dots; RCWA simulation: red curve). The period is 400 nm, the slit width 110 nm 
and the slit height 170 nm. (b) Sketch of light diffraction along the surface of a perfect metal slit array versus the SPP wave at a real metal slit array, 
(c) Sketch of the TEM waveguide mode in a perfect metal slit versus the MIM SPP mode in a real metal slit. 



Similarly, the slit resonance (Fig. 6 (c)) with the ideal metal walls is 
due to the interference of TEM wave with free space momentum; 
while that for real metal (Al) the resonance is due to coupled SPP 
modes inside the metal slits, again having higher momentum than 
free space photons and therefore occurs at longer wavelength. This 
effect can be seen more clearly with Ag metal, which is a better 
plasmonic metal in the visible band. The detailed discussion can be 
found in the Supplementary information online, along with the cal- 
culated dispersion of various cases described above. 

Since the theoretical model is based on the perfect metal approxi- 
mation, all the analysis becomes more accurate at longer wavelengths 
where real metals approach to ideal metals. Therefore, this work 
provides some guidelines in designing the spectrum filters and polar- 
izers not only in the visible-NIR range, but also in the IR and micro- 
wave bands as well. 

In conclusion, the transition from a spectrum filter (resonant 
transmission) to a polarizer (broadband transmission) based on 
the metal slit array structure can be tuned by the period of the array. 
The distinct transmission peak in the former case is due to a coupled 
slit/grating resonance, which broadens with decreasing period as the 
coupled two resonances become separated, and the slit resonance 
which is intrinsically much broader than the grating resonance dom- 
inates the transmission peak. As the slit resonance diminishes due to 
reduced light confinement, a broadband transmission (a polarizer 
behavior) is reached. Another effect of decreasing period is that the 
transmission peak becomes less angle sensitive due to decoupling of 
the grating resonance (angular dependent) from the slit resonance 
(angular independent). Therefore, the polarizer is naturally angle 
insensitive due to the very small period, while a spectrum filter can 
be made angle insensitive by increasing the light confinement in each 
slit to fight against the broadening of the transmission peak with 
decreasing period. 

Methods 

Fabrication of the nano-slit arrays. The nano-slit arrays with smaller periods (from 
200 nm to 500 nm) are created through standard e-beam lithography and reactive 
ion etching. The slit array with a large period (5 um) is created by focused ion beam. 
Both processes were performed on a Aluminum film (170 nm thick) deposited on a 
fused silica substrate. The slit height remains at —170 nm as all slits perforate the Al 
film. The slit width was kept at —110 nm. 

Transmission and reflection measurements. The transmission measurement was 
done by using an invert microscope (Nikon Eclipse TE300) equipped with a 
spectrometer (Ocean Optics HR4000CG-UV-NIR). The illuminating light from a 
halogen light bulb is collimated before impinging on the sample. On the other side, 



the transmitting light is collected by an objective lens with a magnification of 40 X and 
a numerical aperture (N.A.) of 0.75. Since the grating periods (from 200 nm to 500 
nm) we studied are smaller than most part of the wavelength range (from 380 nm to 
880 nm), the objective lens only collects the zeroth order transmission. The reflection 
measurement was performed on the same platform. In this case, a low- magnification 
(4X) objective lens is used to transmit the illuminating light and to collect the 
reflection light. The small N.A. (0.1) of the objective lens ensures little beam 
divergence. The angle-dependent transmission measurement was conducted with a 
small tilting stage which can be fit into the sampling space of the invert microscope. 
The measured angle range (0-10°) is limited by the tilting stage. 

Simulations. The transmission or reflection spectra shown in red lines in Fig. 2 (b) 
and (c) and Fig. 4 (c) and (d) were simulated by RCWA. The field distribution of E x as 
shown in Fig. 5 (a) was simulated by finite element method with the commercial 
software COMSOL Multiphysics. For both methods, the optical constant of Al was 
taken from Ref. 40 and that of fused silica from Ref. 41. 
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